Self-consistent Green's functions calculation of the nucleon mean-free path 
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The extension of Green's functions techniques to the complex energy plane provides access to 
fully dressed quasiparticle properties from a microscopic perspective. Using self-consistent ladder 
self-energies, we find both spectra and lifetimes of such quasiparticles in nuclear matter. With a 
consistent choice of the group velocity, the nucleon mean-free path can be computed. Our results 
indicate that, for energies above 50 MeV at densities close to saturation, a nucleon has a mean- free 
path of 4 to 5 femtometers. 
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The mean-free path. A, of a nucleon in the nuclear 
medium is a basic transport coefficient, essential ingre- 
dient of several theoretical and experimental considera- 
tions. Cross section calculations within Glauber theory 
and transport simulations, for instance, rely directly or 
indirectly on in-medium mean-free paths in an energy 
range from above the Fermi surface to a few hundred 
MeVs [ll-Q. Below the Fermi energy, large values of 
the nucleon mean- free path provide a justification for the 
shell model . At the Fermi surface itself, the existence 
of quasiparticles (qp) is validated by an infinite value of 

Strong interactions between nucleons and the many- 
body correlations they induce play an essential role in 
nuclear systems Their existence precludes the ap- 
plication of mean-field or density functional techniques 
in the calculation of A. As a matter of fact, the fully 
Pauli-blocked nature of these approximations leads to an 
infinite mean-free path below the Fermi surface. More 
sophisticated theoretical approaches have been used in 
the past 0- 12 1 , but even then the hole mean- free path is 



not necessarily well defined [13|, |14|. Moreover, calcula- 



tions have been generally performed at energies arbitrar- 
ily close to the real axis and, as a consequence, ad hoc 
nonlocality corrections have to be introduced, generally 
through the fc— mass approximation Q. 

Many-body Green's functions deal explicitly with 
time propagation of particles and holes in interacting 
fermionic systems, constituting a natural framework to 
compute qp properties [1, 15 1. Pauli principle and 
nucleon-nucleon (NN) correlations are fully taken into 
account. As we shall see in the following, the extension 
of these techniques to the complex energy domain pro- 
vides a more consistent approach to compute qp quanti- 
ties, without further approximations 0, Direct 
calculations on the complex energy plane have been per- 



excitations in solid state applications [l8|. In this Letter, 
we devise such an extension for the case of nuclear matter 
and compute qp spectra and lifetimes, which eventually 
lead to an evaluation of the nucleon mean-free path from 
a fully microscopic standpoint. 

In the past, the lack of microscopic propagators in 
nuclear systems has hampered the calculation of trans- 
port coefficients using Green's function techniques, un- 
like other approaches [l^. The recent implementation 
of the ladder approximation within the self-consistent 
Green's functions (SCGF) framework, however, gives ac- 
cess to the full off-shell energy and momentum depen- 
dence of one-body (IB) propagators. We refer the reader 
to Refs. [2^ [21} for details, but mention here that the 
method takes into account short-range and tensor corre- 
lations. Three-body forces (3BF) are included effectively 
via an average over a third, correlated nucleon |21j| . 

Let us discuss the extension of the Green's functions 
formalism to the complex energy domain, before exam- 
ining the calculation of qp properties. The propagation 
of an excitation in nuclear matter is described by the re- 
tarded propagator, QR{k,t) = 6(<) ({a(fc, t), a'''(A;, 0)}). 
In a uniform system in thermal equilibrium, Qn only de- 
pends on the time difference, t, and the momentum mod- 
ulus, k. Using the Lehmann representation, one finds the 
retarded propagator in energy space: 



gR{k,uj) = 



dcj' A{k,uj') 

llX ijj^ — ijj' 



(1) 



where w+ = cj + i?; with ?; — > 0. All IB operators can be 
built from the spectral function, A{k,ijo). The previous 
expression suggests an extension of the propagator to the 
complex energy plane, replacing a;+ by z: 



Q{k,z) 



dw' A(k,uj') 
27r z — u' 



(2) 



formed in electronic systems since the early 1960's [17 1 



and have been recently employed to describe microscopic 



Such an extension exists, is unique, and gives rise to 
a complex variable function that is analytic off the 
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FIG. 1. Imaginary part of the propagator in the complex 
energy plane, z = E + iF, at zero momentum for a CD-Bonn 
self-energy at p = 0.16 fm'^ and T = 5 MeV. The left plot 
corresponds to the usual propagator, Q{k,z), while the right 
plot represents its analytic continuation, Q{k,z). Solid lines 
show the imaginary part of the propagator just above and 
below the real axis, itA{k,(jj)/2. 



real axis [22j. Close to this axis by above (below), 
Q{k,z) becomes the retarded (advanced) propagator, 
Q{k,uj±) = Qji^ji{k,uj). An illustration of the com- 
plex energy dependence of the IB propagator is given 
in the left plot of Fig. [TJ The discontinuity of Q across 
the real axis is determined by the spectral function, 
lm{g{k,u:^)-g{k,oj+)}^Aik,Lu). 

In the complex plane, the IB propagator fulfills a 
Dyson equation 



gik,z) = 



1 



z — E(fc, z) ' 



(3) 



with a self-energy E which is extended to complex ener- 
gies in analogy with Eq. ^ : 



S(fc,z)= J 



duj 7(fc, u) 

llX Z — ijj 



(4) 



Our strategy will consist in computing ^(k,iJ) from full 
off-shell ladder self-energies and using Eqs. (jl]) and ([3|), 
to extend E and Q to the complex plane. 

Let us now illustrate how the mean-free path can be 
extracted from Q{k,z). One computes A(fc), from the 
nucleon inverse lifetime T{k) and group velocity v{k), 
via A(fc) — v{k)/\T{k)\. Such quantities appear in the 
asymptotic form of the retarded IB propagator [sl. [l5l|: 



t»r- 



->-z7y(A;)e-'^W*e-|rWI* 



(5) 



While r]{k) represents the strength of the excitation, the 
qp spectrum e{k) determines the oscillation frequency of 
the propagator. The spectrum is closely related to the 
group velocity 



v{k) 



de{k) 



(6) 



dk m*{k) 

which is often studied in terms of an effective mass m* (fc) 



It is natural to derive Eq. ([5]) by Fourier transforming 
Qr, under the assumption that Q has a pole of order 1 
in the lower half-plane 1^ 23|. The single pole descrip- 



tion is particularly attractive due to its simplicity and 
its relation to Fermi liquid theory 0]. However, as we 
have seen, Q is analytic for complex z and, as a matter of 
fact, does not have a pole in the complex plane. To find 
Eq. ([5]) , one actually trades the discontinuity of Q across 
the real axis for a pole in the lower half-plane. The pole, 
however, is not associated to Q, but rather to its analyt- 
ical continuation (ac), which is continuous across the 
real axis but nonanalytic in the lower half-plane [2^. Q 
is computed from the complex Dyson equation, Eq. ([3]), 
with E replaced by the ac of the self-energy, E: 



E(fc,2) 



E(fc,z), Im z > 
E*(fc,z), Im z < 



(7) 



a function that is analytic everywhere in the complex 
plane. The right plot in Fig. [T] shows the imaginary part 
of g{k = 0, z) for a CD-Bonn SCGF self-energy, g is 
analytic across the real axis, but develops an isolated 
pole in the lower half-plane. The position of this pole is 
given by the complex equation: 

z{k) = — + ReE(fc, z{k)) + am E(fc, z(k)) . (8) 
2m 

The solution, z{k) = e{k)+ir{k), gives access to the fully 
dressed qp spectrum and inverse lifetime. As we have ac- 
cess to the ac of SCGF self-energies in the complex plane 
via Eq. (jl)) , we are able to compute fully dressed spectra 
and lifetimes for different momenta, densities, and tem- 
peratures. Note that the solution to the previous equa- 
tion need not be unique in the most general case. In our 
nuclear matter calculations, however, we have not found 
any signature of multiple solutions. 

Previous calculations have relied on solving Eq. (jS]) 
using successive approximations for the complex energy 
dependence of E I9l 4lll. At the lowest order, known as 
first renormalization [l8], one completely neglects the de- 
pendence on the imaginary part of z. This provides the 
usual definition of a qp: 



ei{k) = — +Re^{k,ei{k)), (9) 
Zm 

ri(fc) ==Im E(fc,ei(/c)), (10) 

which usually coincides with the peak of the spectral 
function. A second renormalization qp pole is obtained 
by expanding the self-energy around Zi(fc) to first order 
in the imaginary part of z, 



£2(fc) = ei(fc) -Im E(fc,ei(fc))Im 



r2(fc) = ri(fc)Re 



1 



1-E'(zi(fc)) 



1-E'(zi(fc)) 



(11) 



(12) 
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Energy, (n-|i [MeV] Energy, to-ji [MeV] Energy, a>-\i [MeV] 

FIG. 2. Upper panels: spectral function at p — 0.16 fm~^ 
and r = IVIeV for the CD-Bonn interaction. Lower panels: 
absolute value of G in the same conditions. The fully dressed 
pole is indicated by a cross, while the circle (square) show the 
position of the first (second) renormalization quasiparticle. 



In the context of nuclear physics, it has generally been 
assumed that the dependence of S on the imaginary part 
of z is soft and can be ignored in the previous derivatives 
This gives rise to a slightly different qp pole: 



£2'ik) = ei{k) 

r2'(fc) = ri(fc) 



1 -Rel]'(ei(fc)) 



(13) 
(14) 



As we shall see, this approximation is well justified only 
above kp- 

In the following, we present our fully dressed results 
and compare them to previous approximations at p = 
0.16 fm-3. The upper panels of Fig. [2] show the SCGF 
spectral function, as a function of energy, for three differ- 
ent characteristic momenta (k = 0, /cf, and 2^^)- These 
have been obtained from a T = CD-Bonn self-energy 
2l|. The lower panels give the absolute value of the 
analytically continued propagator. Contour levels unam- 
biguously demonstrate the existence of a pole in Q. The 
location of the fully dressed pole is consistent with the 
numerical solution of Eq. ([5]), shown with a cross. Dif- 
ferences between this pole and the first or second renor- 
malization properties are visible at fc = 0. At and above 
the Fermi surface, discrepancies disappear and the fully 
dressed pole coincides with first and second renormaliza- 
tions. This points towards a very soft dependence of S 
on the imaginary part of z for k > kp- Note that, at the 
Fermi surface, calculations yield a zero width, providing 
a verification of Fermi liquid theory from a self-consistent 
perspective 3|. 

Nuclear many-body calculations are subject to uncer- 
tainties associated to the underlying NN interaction as 
well as to the approximation scheme itself. To assess 
them, we summarize in Fig. [3] the results obtained with 
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FIG. 3. quasiparticle properties at p — 0.16 fm~ and T — 5 
MeV for different approximations and two NN potentials: 
CD-Bonn (left) and Avl8 (right panels). Upper panels cor- 
respond to the qp spectrum, central panels to inverse qp life- 
times, and lower panels to effective masses. The different 
approximations are explained in the text. 



two different phase-shift equivalent potentials, the CD- 
Bonn |2j| and the Argonne vis (Avl8) [25*1 interactions, 
ai p = 0.16 fm"'^ and a finite, but rather small, temper- 
ature of r = 5 MeV. The upper panels show the three 
approximations to qp spectra discussed earlier as a func- 
tion of momentum. While above the Fermi surface the 
agreement between all approximations is good, below kp 
the fully dressed pole (solid line) is always more attrac- 
tive than £i (dashed line). In contrast, the second renor- 
malization spectrmn (dash-dotted line) is more repulsive. 
This indicates that successive renormalizations might not 
yield results closer to the fully dressed pole. The inverse 
qp lifetime, shown in the central panels, is bell shaped. 
Remarkably, below the Fermi surface the lifetime is finite. 
Close to A;_F, its absolute value becomes small, but not 
zero due to thermal correlations 0] • Although not shown 
here, we have found that the effect of 3BF at this den- 
sity is small in all the quantities shown [2l|. In contrast, 
many-body approximations other than GF's would yield 
rather different results. Within the Brueckner-Hartree- 
Fock approximation, states below the Fermi surface are 

completely blocked and F = for k<kp [H,!!!. 

To finish the calculation of the nucleon mean-free 
path, a consistent determination of e, F, and v, via 
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upper panel shows the mean-free path obtained with a 
fully realistic self-energy based on the CD-Bonn interac- 



50 100 150 200 250 300 350 
e-H [MeV] 

FIG. 4. Mean-free path of a nucleon in nuclear matter as a 
function of energy. Upper panel: results obtained with a CD- 
Bonn-|-3BF self-energy at T = MeV. The different approx- 
imations are commented on the text. Lower panel: results 
obtained mean-free path from the fully dressed pole for dif- 
ferent NN forces and two different temperatures. The shaded 
band and solid dots correspond to the experimental results of 
Refs. [2, HH], respectively. 



Eq. (|n]), is needed. Let us, for instance, consider the 
nuclear physics renormalization, Eqs. and (ITil) . 

The prefactor on the inverse lifetime is the inverse of 
the oj— mass. The group velocity involves the full ef- 
fective mass, ^ = As a consequence, the 
mean-free path is only renormalized by the fc— mass, 
A2'(A:) = ^Aq, with respect to the uncorrected mean- 
free path, Ao(fc) — /c/[2mlm E(fc,ei(/c)+)]. Thus, con- 
sistency in the spectrum and the lifetime are needed to 
obtain a nonlocality correction . Similarly, the cal- 
culation of the self-consistent mean-free path relies on an 
effective mass computed from the fully dressed spectrum 
rather than on mj. 

The lower panels of Fig. [3] show the effective masses as- 
sociated to the different approximations. Corresponding 
to a more attractive spectrum, the effective mass asso- 
ciated to the fully dressed pole (solid line) is lower than 
that of the usual qp approximation (dashed line). The 
latter is the product of the w— (double-dot dashed) and 
the fc— masses (dotted). Our results confirm the well- 
known behaviors for these quantities: m^^ peaks around 
kp, while TOfe dips in this region. Note, however, that 
for the softer interaction (CD-Bonn), — becomes 1 at 
lower momenta than for the harder force (Avl8). For 
hole states, the effective mass associated with ei shows 
a substantial increase as fc -> 0, associated with the flat- 
tening of the spectrum in that region. 

A summary of final results is presented in Fig. 21 The 



tion supplemented with and Urbana-type 3BF [2l|. As 
expected, we find that the largest differences between 
approximations occur for hole energies, below —20 MeV, 
in a region where it is already relatively small. In con- 
trast, above 50 MeV, all approximations give similar re- 
sults, except for Aq (dashed line), which is not corrected 
for nonlocality and thus should not be taken as a re- 
alistic prediction. A2' (dash-dotted line) is only some- 
what larger than Aq because of the small associated 
with the SCGF results. The kinetic theory prediction, 
A ~ (pcrnp)"^ (dotted lines), is well below all quantum 
in-medium mean- free paths. The latter fiatten at high 
energies, and remain constant, at a value of around 4 — 5 
fm. 

The lower panel of Fig.|4]focuses on the NN interaction 
and temperature dependence of our results. The T = 
mean-free path with 3BF (solid line) is slightly larger 
than that obtained without 3BFs (dashed). The effect 
of temperature is relevant in an area of about 20 MeV 
around the Fermi surface, where the mean-free path is 
finite, although still large. The fully correlated results 
agree with experimental estimates jll, [l^l and suggest 
that A ^ 4 — 5 fm above 50 MeV. The spread between 
different lines is an estimate of theoretical uncertainties, 
which amount to less than 1 fm at those energies. 

To summarize, we have devised a new method to ob- 
tain the mean-free path of a nucleon in the medium. The 
method involves the extension of Green's functions tech- 
niques into the complex plane. The pole of the prop- 
agator gives access to fully dressed qp properties. The 
renormalization induced by this procedure is relevant for 
hole properties. Our approach provides a validation for 
previously used approximations by taking into account 
the full dependence on the imaginary part of the energy. 
With all many-body corrections properly implemented, 
we obtain a mean-free path of around 4 — 5 fm at satu- 
ration density and energies above 50 MeV. Future work 
will systematically assess the density, temperature, and 
isospin asymmetry dependence of the mean-free path. 
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